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PACS numbers: 13.25.Gv, 14.40.Pq 



Two-photon fusion events provide a useful production 
mode to study charmonium states with quantum num- 
bers JP"^ = 0±+, 2±+, 4±+, 3++, 5++, ... 0. 
Dipion transitions among these states have been ex- 
perimentally studied only in one case [4], in contrast 
to the narrower vector states, where dipion transitions 
have been studied extensively. In particular, the transi- 
tion amplitude for ■qd'iS) — ?> ric'n:^'n:~ [3] is expected [5)] 
to have the same approximately linear dependence on 
the invariant-mass-squared of the dipion system as the 
■>Jj(2S) — )> J/T/jTr+TT" decay @. Phase-space integration 
of the squared amplitude, evaluated for the peak masses 
Mjj^ and M^^(^2S) Vi of the rfc and t]c{2S), respectively, 
yields r(r/c(25) ri^'K+Tr-)/r{^j{2S) J/^n+n-) w 
2.9. This leads to the branching fraction prediction 
B{r]ci2S) -> rjcTT+TT-) = (2.2tJ-j)%, where the uncer- 
tainty is due to the uncertainty on the width of the 
r]c(2S) Q. This decay may be further suppressed due 
to the contribution of the chromo-magnetic interaction 



to the decay amplitude Q- 

In recent years, experiments have reported evidence 
for charmonium-like states, such as the X(3872) and 
y(4260) [9], which do not fit well into the conventional 
cc picture. This has prompted much theoretical activ- 
ity and proposals for new models [10]. Several studies of 
these states have been performed with the J/ipn^Tr" fi- 
nal state Uj, but no search using the yycTr+vr^ final state 
has been conducted. Such a search may shed light on 
the quantum numbers or the internal dynamics of these 
states. In particular, it has been suggested [l^l that if the 
X(3872) is the 1^D2 state ?7c2, then the branching frac- 
tion B{X{3872) — > ricTT~^TT~) could be significantly larger 
than i3(X(3872) — > J/^tt+tt^). The quantum numbers 
jP'-^ = 2 ^ of the 77c2 are consistent with the results 
of an angular analysis of X(3872) — >■ J/i/jTr+Tr" fl3'] and 
would allow production of X (3872) in two-photon fusion. 

We present herein a study of the process 77 — > X — ^ 
77c7r+7r~, where X is one of the resonances Xc2{^P), 
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7]c(2S'), X(3872), X(3915), or Xc2(2P), and the 77^ is re- 
constructed in the final state K^K'^'it~ 

The data sample was collected with the BaBar detector 
at the PEP-II asymmetric-energy e"''e~ collider located 
at the SLAC National Accelerator Laboratory. It consists 
of 429.1 ± 1.9 fb~^ collected at the energy of the T(45') 
resonance, constituting the entire BABAR T(45) dataset, 
and 44.8 ± 0.2 fb"^ collected about 40 MeV below the 
T{AS) resonance. The BABAR detector is described in 
detail in Ref. [11]. 

Samples of Monte Carlo (MC) simulated events are 
analyzed with the same reconstruction and analysis pro- 
cedures as the data sample, following a GEANT4-based [13] 
detector simulation [l^]. Simulated background samples 
include e+e~ qq continuum events (g = u, d, s, c) gen- 
erated with JETSET [18*], T{4S) -> BB decays gener- 
ated with EvtGen pJi] and JETSET, and e+e" ->■ t+t' 
events generated with KK . In order to study initial- 
state-radiation (ISR) background and the invariant-mass 
resolution, a sample of e^e~ — >■ ^ip(2S) events with 
ip{2S) J/?/)7r+7r^ and J/ip — >■ K^K~^tt~ is generated 
with EvtGen. The GAMGAM ^ generator is used to 
generate signal event samples for each of the X states 
studied, with the decay X — > tIcTt^t:^ simulated with an 
amplitude that is uniform throughout the decay phase 
space, independent of the final-state kinematic variables. 
The decay tJc — > K'^K^tt" is generated with a uniform 
amplitude or with equal and incoherent i^Q (1430)~i^+ 
and Ko*(1430)°X° contributions. The GAMGAM gen- 
erator is also used to generate 77 — > 77c — > K'^K^tt^ 
events. 

The analysis is performed with two data samples. The 
sample used to search for the process ^ X ~^ rjcn^n^ 
is referred to as the "main sample" . Properties of the r/c 
and its decay into K^K^Tr~ are studied with a separate 
"control sample" of ^ rjc ^ K^K^n^ events. For 
the main (control) sample, we select events that contain 
six (four) charged-particle tracks. 

For both samples, charged kaon candidates are iden- 
tified using likelihood values calculated from measure- 
ments of specific energy loss and information from a de- 
tector of internally reflected Cherenkov radiation. All 
other tracks are assumed to be pions. A candi- 
date is reconstructed by fitting a tt+tt" pair to a com- 
mon vertex, with invariant mass in the range 0.491 < 
m{'K^n~) < 0.503 GeV/c^. A kinematic fit is performed, 
constraining TO(7r+7r~) to the nominal ifg mass ^22i] . An 
77c — >■ K^K^TT~ decay candidate is reconstructed by com- 
bining a candidate with a and a it~ and re- 
quiring the resulting invariant mass to lie in the range 
2.77 < ■m{KlK+'K-) < 3.22 GeV/c^. In the main sam- 
ple, the decay X — 77c7r"'"7r~ is reconstructed by com- 
bining an r]c candidate with the remaining two tracks in 
the event. A kinematic fit is applied, requiring the X- 
candidate decay vertex to be consistent with the e^e^ 
interaction region. The angle aj^o between the mo- 



mentum vector and the line connecting the 77c and the 
Kg decay vertices is required to satisfy cosa^o > 0.99. 

In the control sample, we require the polar angle (the 
angle with respect to the beam axis) 0,,^ of the rjc can- 
didate to satisfy |cos6',,^| > 0.99 and the transverse mo- 
mentum of the 77c candidate to satisfy < 0.5 GeV/c, 
both in the center-of-mass (CM) frame of the e+e^ sys- 
tem. The extra energy in the event, defined as the to- 
tal energy in calorimeter clusters not associated with the 
identified tracks, is required to satisfy i?cx < 0.5 GeV 
in the CM frame. The m{K'^K~^7r~) distribution of the 
selected control-sample events, shown in Fig. [Tfa), ex- 
hibits clear 7;c and J/ijj peaks, with the J/ip produced in 
ISR events. In the main sample, continuum background 
is strongly suppressed with the requirements ]cos0xl > 
0.85, < 1.5 GeV/c, and Sox < 0.8 GeV, where cosOx 
and are the polar angle and transverse momentum of 
the X candidate. In addition, the total visible energy in 
the event, obtained from all charged tracks and calorime- 
ter clusters, is required to satisfy Eyis < 10 GeV in the 
laboratory frame. 

In the main sample, we suppress QED background by 
requiring R2 < 0.7, where R2 is the ratio of the sec- 
ond and zeroth Fox- Wolfram moments 23|. We sup- 
press background due to ISR events with a requirement 
on the missing mass squared rn^-^^^ = (jPe+e- ~ VxY > 
10 GeV^/c^, where Pe+e- {Px) is the total 4-momentuni 
of the beam particles {X candidate). 

Additional background suppression in the main sample 
is obtained by using the Dalitz plot for the t^c candidates. 
The Dalitz plot is shown in Fig. Hfb) for control-sample 
events in the t^c peak region 2.94 < m{KlK^TT^) < 
3.02 GeV/c^, and in Figs. \V[c) and [Hd) for main- 
sample events in the lower and upper m{K'^K~^'!T~) side- 
bands 2.8 < m{K'^K+TT-) < 2.9 GeV/c^ and 3.05 < 
m{K'^K+TT~) < 3.2 GeV/c^, respectively. These dis- 
tributions indicate that true 77c — K';iK^7r~ decays 
often proceed via intermediate i4rQ(1430) states, while 
background events contain K*{892) decays and ran- 
dom combinations. Taking advantage of this differ- 
ence to suppress non-7yc background in the main sam- 
ple, we require \m'^{K^'K~) ~ Af|..(^43o)_ ] < 0.5 GeV^/c"* 
or Im^iK+TT-) - M|..(^43o)ol <''o.5 GeVVc*, and ex- 
clude events that satisfy \m'^{K'^7r^) - 
0.35 GeVVc^ or \m^{K+Tr-)-M]^,^^g^^a 
where Af^ is the peak mass of resonance R [2^. The 
Dalitz-plot region selected by these criteria is enclosed 
within the solid lines in Figs, [ijb), [T{c) and[ljd). The 



^■^^•(892)- I < 

< 0.2 GeV7c^ 



criteria are the result of maximizing / \J ^ where 

= (63.5 ± 3.2)% is the efficiency of the Dalitz- 
plot requirements for 77c decays, determined by fitting 
the m{K^K^'K^) distribution of the control-sample, and 
^SB ~ (30.74 ± 0.21)% is the corresponding efficiency 
for main-sample events in the m{K^K^T:^) sidebands. 
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The -D° — >■ ii'+TT" band, evident in Fig. [U^c), becomes 
insignificant following a neural- network requirement, de- 
scribed below. Therefore, no explicit effort is made to 
remove this source of background. 



■^200 

u1000 

C5 

uj 800 

o 

9 600 
o 

« 400 

§i 200 
UJ 







(a) 










d) 5- 








O 








E 3r 








2r 








1r 



2.8 



2.9 3.0 3.1 3.2 

m(K''Kjt) (GeV/c^) 



(b) 



> 
(5 

H 

E 







(c) 






4H 
















n 






o 



7 

0) 6 
t5 

E 4 



3S 



(d) 



3 4 5 6 

m2„ (GeV^/c") 



% 1 2 3 4 5 6 7 8 



m^,^ (GeV^/c") 



in 0.5 
o 
c> 

^0.4 

0) 

■>-0.3 
> 

■■S 0.2 



I I Signal 

Background 



(e) 





2.8 2.9 3.0 3.1 3.2 

mCK^Kii) (GeV/c^) 



FIG. 1: (a) The m{KsK+n~) distribution for the control 
sample. The vertical lines indicate the rjc peak mass region. 
Also shown are the KsK'^tt~ Dalitz-plots for (b) control- 
sample events in the rjc peak mass region and for main-sample 
events in the (c) lower and (d) upper rjc mass sidebands. Solid 
black lines indicate the regions defined by the Dalitz-plot se- 
lection criteria. The dotted blue box in the upper left corner 
of (c) and (d) indicates the Dalitz-plot-sideband background 
region used for the neural-network training, (e) The neural- 
network output- variable distributions for the Dalitz-plot side- 
band (hatched) and signal MC. (f) The result of the step-1 
fit (see text). 

Further background suppression is achieved by com- 
bining six variables into a neural-network discriminator. 
Two of the variables are Ecx and pj^. The other four 
variables, each of which can take on one of five discrete 
values, are the outcomes of kaon- and pion-identification 
algorithms applied to each of the four charged-particle 
tracks that are not the daughters of the K'^ candi- 
date. The neural network is trained with main-sample 



background events in the Dalitz-plot sideband region 
m^{K+Tr-) > 2.5 GeVVc^, m^K^n-) < 1.5 GeV7c^ 
indicated by the dashed boxes in Figs. [TJc) and [Hd). 
This region is chosen since it contains only (3.40 ±0.66)% 
of r]c decays in the control sample. We find only in- 
significant differences in the neural network signal-to- 
background separation when using different signal sam- 
ples or the mirror Dalitz-plot region m'^{K^TT~) < 
1.5 GeVVc", m^{K'^Tr~) > 2.5 GeV^c* for the back- 
groimd. 

The distributions of the output-variable Vnn are shown 
in Fig.[TJe). We find the optimal selection on this variable 
to be T^^N > 0.84. The efficiency of this selection is 72% 
for the rye (25'), and varies by up to 4%, depending on the 
X mass. The background efficiency is (10.4 ± 0.2)% for 
the neural-network training region and (7.4 ± 0.2)% for 
the mirror region. 

We find 2863 main-sample events that satisfy all the 
selection criteria, with only about 700 events expected 
from non-77 background MC. We conclude that the ma- 
jority of the background is due to 77 events, for which we 
have no generic generator. More than one X candidate 
is reconstructed in 3.8% of the events. In these cases, we 
select the candidate for which m'^{K'^Tr^) or m^(i4r+7r~) 
is closest to the i<rQ(1430) peak. 

In addition to these samples, an ISR- produced sam- 
ple of ip{2S) J/f/jTr+Tr" events is used to evaluate a 
systematic uncertainty associated with the detector res- 
olution. This sample is selected in the same way as the 
main sample, except that the neural-network and Dalitz- 
plot selections are not applied, the KgK'^Tr~ invariant 
mass is required to be between 3.0 and 3.2 GeV/c^, and 
m^jgg must be less than 1 GeV^c**. 

We define four categories of events in the main sample: 
signal corresponds to 77 — > X — > TycTr+Tr" events; com- 
binatorial background (CB), which is by far the most 
copious background, arises from random combinations of 
final-state particles; events with a true rjc K'^K~^7r~ 
decay and two pions not originating from an X reso- 
nance decay are categorized as r^c-peaking background 
{ricB); X-peaking background (XB) corresponds to de- 
cays X K'^K'^TT^TT^TT^ that do not proceed through 
an intermediate 7]^. 

The extraction of the signal yield proceeds in two steps. 
In step 1, we determine the values of the m{KgK~^TT~)- 
distribution parameters of the combinatorial background 
from a one-dimensional fit to m{K'^K'^ Tr~), without any 
restrictions on ■m{K'^K'^7T^TT^TT^). 

In step 2, we extract the signal yield for each X 
resonance hypothesis from a two-dimensional fit to the 
■m{K'^K^TT^) versus m{K'^K^Tr^TT^TT^) distribution for 
events in an m{KgK~^TT~Tr^TT^) window around the res- 
onance peak. The fits use the unbinned, extended- 
maximum-likelihood method and are performed with the 
RooFit package [iij . 

From events in the 'm{K'^K^Tr^) sidebands, we ob- 
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serve that all correlation between the ■m{K'^K~^TT~TT^TT~) 
and 'm{K'^K^ Tr~) distributions for the combinato- 
rial background is accounted for by the phase space 
<i>{Tn{K°K+Tr-),m{K°K+Tr-Tr+Tr-)) of the three-body 
final state consisting of the 7r+, tt^, and the {K'^K'^tt") 
system. This is used to construct the probability-density 
function (PDF) ^ of the step-1 fit. This PDF is a 
function of ■m{K'^K~^TT~) with m{K'^K~^TT~TT^TT~) as a 
conditional variable, and is given by: 

'^(malms) = Nr^^n^^im^) + iV^^H^j^ (malms), (1) 

where N^^ i^iia) t'^^ number of events with (without) 
a true rjc — J> K'^K^tt' decay. We have used the notation 
ma = rn{K'^K^TT^) and = ■m{K'^K^Tr^Tr^Tr^) for 
brevity. 

The PDF for non-yyc events in Eq. ([T]) is 

H,^^ (malms) = P2(m3; ai, 02, m3)$(m3, ms), (2) 

where V2im3; ai, 02, m!]) — 1 + ai(m3 — m!]) + a2(m3 — 
ma)^ is a second-order polynomial and m!] = 3.0 GeV/c^. 
Determination of the coefficients ai, 02 is the main pur- 
pose of the step-1 fit. The PDF for rjc events in Eq. ([T|) is 
^»)c(™3) — W(m3; Af^^, r^^, rvig), where W is a relativis- 
tic Breit-Wigner function [{ml - M^J^ + M^T^^]'^ 
convolved with a detector resolution function TZ{m3 — 
m3;7Vi3) that depends on a set of parameters r,„3 and 
the true invariant mass 7713 of the K^K~^tt~ system. 
The resolution function is the sum of two Crystal Ball 
functions [2^ with oppositely-directed tails and common 
Gaussian-parameter values. The resolution-function pa- 
rameter values are determined from a fit to the MC. 

In addition to ai and 02, the parameter values deter- 
mined in the step-1 fit are the yields iV^^ and iV^^, and 
the mass M^^ and width F^^ of the rjc peak. In order 
to obtain M^^ and F,,^ from the data, the step-1 fit is 
performed simultaneously for the main sample and the 
control sample. The PDF for the control sample is 

H'ims) = N'jf^,W{m^;Mj/^,Tj^,r^^) (3) 

Additional control-sample parameter values determined 
in the fit are the peak J ftp mass Mj/^, the background 
parameters a'l, Oj, and the rjc, J/V', and background 
event yields N'^^, N'j^, and iV^gj. 

The m{K2;K'^TT~) distribution of the data and the 
step-1 PDF are shown in Fig.[Tl(f). The fitted parameter 
values are ai = 1.24 ±0.19 (GeV/c^)"!, 02 = 0.2 ± 1.4 
(GeV/c2)-^ 7V,,^ = 50 ±37, N'^^ = 10350 ± 300, and 
^J/ip ~ 1877 ± 90. The large relative uncertainties 
for ai and 02 are the result of the near linearity of 
the m{KgK'^7T~) distribution and the correlation be- 
tween the two parameters, which is taken into account 
in the evaluation of systematic uncertainties. The rjc pa- 
rameter values determined in the step-1 fit are F^^ = 



31.7 ±1.5 MeV/c2 and M^^ = 2.98285 ±0.00038 GeV/c^, 
where the uncertainties are statistical only. These results 
are consistent with previous measurements [22j . 

The PDF for the step-2 fit is a linear combination of 
the PDFs of the four event types, 

r = A^sig^ig + iVcB^CB + iV„.B^r,.B + NxB^XB- (4) 

The signal PDF is a relativistic Breit-Wigner function 
convolved with the resolution function, for both m3 and 
ms: 

^sig(m3,m5) ^H,^^{m3)W{m5;Mx,Tx,r„i^), (5) 

where Mx and Tx are the known mass and width of 
the resonance of interest [ll, [l^ , and fm^ are the pa- 
rameters of the m(Ar°i(r+7r~7r"'"7r~) resolution function 
7?.(ms — ms; rm^), obtained from a fit to signal MC. The 
combinatorial-background PDF is 

7'cB(m3,ms) -%(m3|ms)C2(ms;6f^,6^^), (6) 

where C2(ms; &f ^, &2 ^) is a second-order Chebychev 
polynomial with first- (second-) order coefficients 
(&2^). The 77c-peaking background PDF is 

^>)eB(m3,ms) ^Hr,,{m3)Ci{m5;bf^), (7) 

where Ci (ms ; bi"^) is a first-order Chebychev polynomial. 
The A-peaking background PDF is 

VxB{m3,m5) = 7'i(m3;cfB,m°) (8) 

where Pi(m3; c;^^, mjj) is a first-order polynomial. The 
parameter values determined with the step-2 fit are the 
four yields of Eq. (jH) and the background shape param- 
eters 6?^, 6^3, 6?=^, and cf 

The step-2 fit is performed four times in different 
m{K^K~^Tr^n^Tr^) windows, fitting for the (1) Xc2{^P), 

(2) r]c{2S), (3) X(3872) and A:(3915), or (4) A:(3872) 
and Xc2(2P) resonances. A simultaneous fit to the three 
resonances A (3872), A (3915), and Xc2{2P) is observed 
to be unstable when tested with parametrized MC exper- 
iments, due to the large number of fit parameters, small 
signal, and large overlap of the A(3915) and Xc2{2P) 
lineshapes. Therefore, we conduct fits (3) and (4) sepa- 
rately to test for the existence of a signal for either set 
of lineshape parameters. The m{KgK'^Tr^'K^Tr^) and 
m(A'gA'+7r~) distributions and fit functions are shown 
in Fig. O The difference between the fit function of fit 

(3) and that of fit (4) is almost indistinguishable within 
the thickness of the curve in Fig. [HJf). The fitted signal 
yields are summarized in Table [H 

No significant signal or peaking background is ob- 
served in any of the fits. However, a hint of A-peaking 
background is visible in the Xc2{^P) and ric{2S) fits of 
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Figs.[2tb) and (d), with event yields of 33±14 and 47±24, 
respectively, where the uncertainties are statistical only. 
This may be due to decays of Xc2{^P) and ric{2S) into 
KgK'^TT~T:'^Tr~ [23|, which are suppressed in this analy- 
sis by the 2.77 < m{K°K+7T-) < 3.22 GeV/c^ require- 
ment. Fits (3) and (4) yield insignificant X-peaking back- 
ground, roughly canceling the negative signal yields. The 
results shown in Table H] for the X{3872) are obtained 
from fit (4). The X(3872) yield from the X(3915) fit 
is 1.6 events lower. Since no signal is observed for the 
X(3915) or the Xc2(2P), we obtain a conservative upper 
limit on the yield of the X(3915) (xc2(2P)) by fixing the 
Xc2(2P) (X(3915)) yield to zero. 
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FIG. 2: Distributions of (a,c,e) m{K°sK+n^) and (b,d,f) 
m{KsK+-K^Tr+TT') with the step-2 fit PDF overlaid for the fit 
regions of the (a,b) Xc2(lP), (c,d) ric{2S), and (e,f) X(3872), 
X(3915) and Xc2{2P). The vertical dashed lines in (f) indi- 
cate the peak mass positions of the X(3872), X(3915), and 

xc2(2P) m- 

We estimate systematic uncertainties on the signal 
yields associated with the fit procedure by repeating the 
fits with the variations described below and adding the 
different uncertainties in quadrature. We account for un- 



certainties in the X mass and width values by varying 
them within their uncertainties |22|. This is the source 



of the largest signal-yield systematic uncertainty, except 
for the Xc2(l^')- The order of the polynomial in each 
PDF is varied to account for uncertainties due to back- 
ground modeling. We vary the resolution-function Gaus- 
sian width by 2 MeV/c^ for the m{K°K+TT-) PDF to 
account for a difference between the J/ip width in MC 
and in the control sample, and by 0.9 MeV/c^ for the 
m{K'^K+n-TT+n-) PDF due to a difference in the ip(2S) 
width between MC and data. An additional uncertainty 
is evaluated by using the sum of three Gaussians to de- 
fine the resolution function. To address the possibility 
that correlations between the and 7713 distributions 
are not taken fully into account by the phase-space fac- 
tor $(7713,7715) in Eq. we replace the parameters 
of Eq. ([2]) by ai{l + a'/m^). The values of the param- 
eters a" are found to be consistent with zero, and we 
conservatively use their uncertainties to evaluate the sys- 
tematic uncertainty on the signal yield. The effect of 
not accounting for phase-space correlations between TO3 
and 7715 in the signal and ?7c-peaking background PDFs 
is determined to be small compared to other systematic 
uncertainties, except for the Xc2(l^'), for which this un- 
certainty is dominant and equals 2.4 events. Statistical 
uncertainties from the step-1 fit are propagated to the 
step-2 fit, accounting for correlations among the param- 
eters. 

We test the entire fit procedure using parameterized 
MC experiments generated with the PDFs of Eqs. ([1]) 
and (U]). A bias of up to two events on the signal yield 
is found and used as a correction that is accounted for 
in the values shown in Table HI A systematic uncertainty 
on this correction is evaluated by repeating this study 
after varying the generated signal yield by its statistical 
uncertainty in the data fit. 

Since the background is mostly combinatorial, we do 
not expect significant interference between signal and 
background. In addition, the small signal yields make the 
evaluation of such interference effects unreliable. There- 
fore, we do not attempt to account for possible interfer- 
ence. 

We evaluate systematic uncertainties on track and K'^ 
reconstruction efficiencies, accounting for the momentum 
and angular distribution of signal tracks, as well as on 
the uncertainty of the Dalitz-plot requirement efficiency. 
A 2% systematic uncertainty is assigned due to differ- 
ences between the distributions of the selection variables 
in the control sample and in 77 — ^ 7;c — K'^K^tt^ MC. 
Differences between the data and MC distributions of the 
particle-identification variables are studied using a high- 



purity sample of D 



*+ 



7r+D°, £)° K-TT+ events, 



and found to have negligible impact on the efficiency. 

We account for uncertainties in the X — s> tIcTt^tt^ 
decay amplitude, which is uniform in our simu- 
lated signal samples, by weighting events according to 
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(m^(7r7r) — 4Af^) Q, where m'^{mT) is the squared di- 
pion mass and is the 7r~ mass. From the weighted 
sample, we extract an efficiency correction of up to 4.6% 
(incorporated into the values in TableHl and a systematic 
uncertainty of the same magnitude. 

Finally, we account for a 0.45% uncertainty on the 
integrated luminosity, for the uncertainties on the ifg, 
77c, and ric{2S) branching fractions 0, and for MC- 
statistical uncertainties. 

The results are summarized in Tabic HI From the signal 
yield iVgig of each resonance, the integrated luminosity 
L, and the signal efficiency e, we compute the product 



iV,ig/(Le) of the e^ 



— > Xe~^e production cross 
?7c7r+7r~ branching fraction. We 



section and the X 
also evaluate the results in terms of the product T^-yB, 
where F^^ is the two-photon width of the resonance, by 
utilizing the GAMGAM generator to determine the cross 
section as a function of F^-^. A 10% uncertainty is as- 
signed to the GAMGAM calculation [1]. Since we find 
no significant signal for the X resonances, we calculate 
90% confidence-level (CL) Bayesian upper limits on these 
quantities, assuming a Gaussian likelihood incorporating 
statistical and systematic uncertainties. 

Using the efficiency-corrected yields for the Xc2 and 
77c (25) from [H and the branching fractions of their de- 
cays into KgK~^Tr~ , we find the relative branching frac- 
tions 

BM2S) 



77c7r ' TT ) 

B{r]c{2S) KOK+n-) 



= 4.9 



+3.5 



3.3 ±1.3 ±0.8, (9) 
= 14.5+oV±7.3±2.5, 



where the first uncertainty is statistical, the second is 
systematic, and the third is due to the uncertainty on 
B{r]c K°K+Tr-) (13 • The 90% CL upper limits on 
the two ratios in Eqs. ([9]) are 10.0 and 32.9, respec- 
tively. Using i3(?7c(2S') -j. K°K+n-) and S(xc2(lP) ^ 
K°K+TT-) from Ref. |22|, we obtain the 90% CL upper 
limits S(77c(25) -)> 7?c7r+7r-) < 7.4% and S(xc2(lP) ^ 
77c7r+7r-) < 2.2%. 

In summary, we report a study of the process 
77 —5- 77c7r+7r^ and provide, for the first time, upper 
limits on the branching fractions of Xc2(l^') and ric{2S) 
decays to ric'K~^iT~ relative to the branching fractions of 
the decays into K^K^'k~ . We also report upper limits 
on the products aB and T^^B for the Xc2(l-P), vA^S), 
X{3872), X(3915), and Xc2(2P) resonances. 
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